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We present incoherent quasi-elastic neutron scattering measurements in a wavevector transler 
range from 0.4 to l.GA~^ on liquid n-hexane confined in cylindrical, parallel-aligned nanochan- 
nels of 6 nm mean diameter and 260 /xm length in monolithic, mesoporous silicon. They are com- 
plemented with, and compared to, measurements on the bulk system in a temperature range from 
50 K to 250 K. The time-of-flight spectra of the bulk liquid can be modeled by microscopic trans- 
lational as well as fast localized rotational, thermally-excited, stochastic motions of the molecules. 
In the nano-confined state of the liquid, which was prepared by vapor condensation, we find two 
molecular populations with distinct dynamics, a fraction which is immobile on the time scale of 1 ps 
to 100 ps probed in our experiments and a second component with a self-diffusion dynamics slightly 
slower than observed for the bulk liquid. No hints of an anisotropy of the translational diffusion 
with regard to the orientation of the channels' long axes have been found. The immobile fraction 
amounts to about 5% at 250 K, gradually increases upon cooling and exhibits an abrupt increase 
at 160 K (20 K below bulk crystallization), which indicates pore freezing. 



I. INTRODUCTION 

The stochastic, thermally-excited motions within 
molecular hquids spatially confined on the nanometer 
scale is still attracting_great interest both in fundamental 
and applied sciences.fj Q In particular, self-diffusion of 
molecular hquids confined in mesoporous hosts plays a 
dominant role in catalysis and adsorptive separation. Q 
It has also been discussed with regard to drug delivery 
applications and the origin of adsorption-desorption 
hysteresis found in mesoporous materials. Q In terms of 
fundamental science the Brownian motion of molecules 
in pores a few nanometers across or in microchannels is 
interesting, since concepts concerning random motions in 
highly confined, crowded geometries @ or for molecules 
exposed to external fields (i.e. interaction potentials with 
the confining walls) can be scrutinized. More- 
over, the direct relation between stochastic motions and 
viscous properties of a liquid (expressed by the Stokes- 
Einstein equation) allows one to explore the fluidity of 
liquids in nano-scale structures and in the proximity of 
the confining solid walls by studies of their self-diffusion 
dynamics. This is of obvious importance in the emerging 
field of nano- and microfluidics, where the exploration of 
these properties can be experimentally very demanding. 

[iiil 

There is a sizeable number of experimental techniques, 
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most prominently nuclear magnetic resonance, gravimet- 
ric uptake measurements, x-ray and light correlation 
spectroscopy as well as quasi-elastic neutron scattering, 
which allow one to study diffusion dynamics in nano- 
and mesoporous matrices. The advent of monolithic 
mesoporous hosts with tailorable channels a few nanome- 
ters across, e.g. mesoporous alumina and silicon mem- 
branes with parallel-aligned channels, have additionally 
increased the interest in this field and extended the ana- 
lytical opportunities, for it allows a better and /or simpler 
comparison of theory and experiment, [l^, [13, [25| 
Many studies of diffusion in restricted geometries are 
aimed at a comparison to the spatially unconfined, i.e. 
bulk, system. Over the years it has turned out, how- 
ever, that this comparison is often hampered, if one has 
to solely rely on diffusion data reported in the litera- 
ture and extracted with other measurement techniques. 
Even for identical bulk systems different experimental 
methods can deliver deviating dynamical properties and 
upon confinement these differences in the measured quan- 
tities can even be more pronounced. (26l - [28j Different 
methods probe different length scales and thus also dif- 
ferent time scales of diffusion processes. Translational 
self-diffusion in a crowded melt is, however, a highly 
cooperative phenomenon and can necessitate large scale 
molecular rearrangements depending on the investigated 
diffusion length. This translates to a dependency of the 
diffusion dynamics on the diffusion length, and thus dif- 
fusion time investigated. This was particularly convinc- 
ingly demonstrated by seminal Molecular Dynamics sim- 
ulations on unconfined, molecular liquids [2^ [s^] and is 
all the more of importance for spatially nano-confined liq- 
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uids, where the bare geometrical restriction can hamper 
both the movements of the diffusing molecule as well as 
the necessary mesoscopic rearrangement of the surround- 
ing molecules. [Sll, [32| 

This study tries to circumvent the challenge with respect 
to the comparability of self-diffusion properties gained by 
different methods by an investigation of both the bulk 
and spatially confined dynamics with an identical tech- 
nique under similar experimental conditions. We report 
Incoherent Quasi-Elastic Neutron Scattering (IQENS) 
measurements on liquid n-hexane in tubular channels 
of 6 nm, that is nine times the length of the molecule. 
Hence, no single-file diffusion behavior is expected, how- 
ever the interaction with the attractive confinin g w alls 
is expected to affect the Brownian dynamics. |33l435l | 
Previous quasi-elastic neutron scattering measurements 
on liquid hexane in mesoporous silica were interpreted 
in terms of pronounced reductions in the translational 
self-diffusion of the molecules in the pore center (by an 
order of magnitude in the translational diffusion coeffi- 
cient) and an even more slowed-down component in the 
pore wall vicinity. [36j This study was performed with 
sponge-like, porous silica gel glasses and hence no inves- 
tigation of anistropic diffusion with regard to the long 
pore axes was possible, an additional goal of this study. 
Liquid hexane was a perfect candidate to study molecu- 
lar motions in a bulk liquid (BL) and spatially confined 
liquid (CL) with IQENS techniques at moderate temper- 
atures. It is a volatile room temperature liquid and in a 
neutron beam its molecules cause mainly incoherent scat- 
tering due to the high content of hydrogen atoms. Other 
scattering contributions from hydrogen or carbon atoms 
are not significant. Porous silicon was the preferred host 
to investigate the CL. It offered a natural environment to 
study the effect of nano-confinement and reduced dimen- 
sionality on stochastic motions. Its tubular/cylindrical 
pores have diameters of only a few nanometers but ex- 
tend to lengths of a few hundred microns therefore creat- 
ing an effective confinement in two dimensions for guest 
molecules like hexane. Its pores are parallel-aligned along 
the < 100 >-direction in single crystalline wafer which in 
combination with the huge length to diameter ratio not 
only raises the question whether molecular motions along 
the pore axis are different from radial motions but also 
offers the unique opportunity to study these potential 
anisotropics in selected scattering geometries as discussed 
in more detail below. 



II. SAMPLE PREPARATION 

Liquid hexane (99.99%) was stored in sealed bottles 
in a dry atmosphere of a glove box. Bottles were only 
removed and opened to prepare BL and CL for the 
experiments. In this way contamination of hexane with 
water, which slowly builds up in a humid environment 
was kept at a minimum. 

Electrochemical anodization in hydrofluoric acid was 



applied to etch porous structures in 500 /im thick single 
crystalline silicon < 100 >-plates. Anodization parame- 
ters as wafer doping, electrolyte concentration, etching 
time and current were adjusted to facilitate the growth 
of 260 fim thick porous layers whith a total porosity 
of 45% at a average pore diameter of 6 nm. Etched 
substrates were stored in a container with distilled water 
to prevent contamination with hydrocarbons and only 
removed for experiments. 

The BL as well as CL sample were prepared directly 
at the beamline. The setup to prepare BL included a 
traditional helium cryostate, a standard sample stick 
and a flat aluminum cell. The indium sealed aluminum 
cell with a width of 30 mm, a height of 50 mm and a 
depth of 0.5 mm contained next to the BL several sheets 
of thin aluminum to reduce the effective thickness of the 
hexane sample to avoid multiple scattering. 
Setup and procedure to prepare the CL sample were 
more sophisticated in order to avoid the formation of 
excess bulk liquid outside the pores in the filling process. 
We employed a sample stick equipped with a stainless 
steal capillary, an air-tight network of small pipes, gas 
containers, pressure gages and valves, subsequently ref- 
ered to as "gas handling" and a high vacuum pumping 
station. The sample cell was tailored to match the 
dimensions of the 30 mm wide, 50 mm high and 0.5 mm 
thick porous matrix. The stainless steal capillary at the 
stick acted as narrow pipeline between the sample in 
the aluminum cell and the external gas handling. The 
gas handling itself served as reservoir and distribution 
system for hexane vapor. 

The sample cell was evacuated at 373 K to 5x 10~^ mbar 
before the pores were stepwise filled with liquid hexane. 
This removed residual water from the inherently hy- 
drophobic porous silicon host. For filling purposes, the 
sample was kept at a temperature of Tq = 273 K, the gas 
handling remained slightly warmer at room temperature 
close to T = 300 K. Then volumetrically controlled 
sorption steps allowed to successively physisorb well 
defined amounts of liquid hexane inside the pore space 
up to complete filling with iVo hexane molecules at the 
saturation vapor pressure of p = Po(Tq) = 61 mbar. 
A complete adsorption-desorption cycle for hexane in 
porous silicon is shown in Fig. [T] It exhibits on ad- and 
desorption branch the number of physisorbed molecules 
N as function of the reduced pressure p/po and is 
commonly referred to as sorption isotherm. Two distinct 
sorption regimes on each branch are due to different 
geometric distributions of hexane molecules in the pores 
depending on the relative pressure. At low pressures, 
that is p/po < Pcad for adsorption and p/po < Pcdc for 
desorption, the molecules physisorb preferentially close 
to the pore wall.[33| In the hysteretic regime above the 
transition pressures Pcad and Pcdo, which relate directly 
to the pore diameter, capillary condensation occurs 
in the pore centre [S^. A detailed discussion of the 
hysteresis (pcad > Pcdc) can be found in a multitude of 
references IS&Mll. 
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FIG. 1: (Color online) Sorption isotherm for hexane con- 
densed in porous silicon. The number of adsorbed molecules 
A'^ is shown as function of the reduced pressure p/po- Open 
symbols refer to adsorption, closed symbols refer to desorp- 
tion. Insets illustrate different geometric distributions of con- 
densate in the pores as discussed in the text. 



For scattering experiments on CL's, the porous host 
was filled at Tq = 273 K up to only 95%, decoupled 
from the gas handling and eventually brought on 
temperature T < Tq. Consequently all experiments 
were performed under quasi-isosteric conditions, and 
scattering contributions from coexisting BL could be 
ruled out definitely. If the pore volume was filled up 
to 95% only at 273 K and the cell was decoupled from 
the reservoir prior to cooling, the filling fraction of 
the host with hexane could not change by more than 
2% through condensation from coexisting vapor while 
cooling. The excess volume of cell and capillary, that 
is the combined volume of cell and capillary minus the 
volume of the empty sample, were neglectable com- 
pared to the pore volume and did not contain enough 
hexane vapor to increase the filling of the host signif- 
icantly. Conditions for bulk formation in the sample 
cell, namely p/po = 1 and N/Nq > 1, were never fulfilled. 



III. SCATTERING EXPERIMENTS 

Scattering experiments were performed at the time- 
of-flight spectrometer IN5 at Institut Laue-Langevin 
(Grenoble). The incoherent quasi elastic scattering of 
cold neutrons {Ei = 2.08 meV, Aj = 0.627 nm) served 
as probe for stochastic motions of hexane molecules in 
the liquid phase. The scattering signal S{Q,uj), that 
is the number of neutrons scattered per second into a 
small solid angle dfl with an energy transfer between 
hu! and h{ijj + d(jj), was recorded as a function of wave 
vector transfer Q and energy transfer fiuj. Wave vector 



transfers Q between 0.4 and 1.6 were acces- 
sible with a spatially fixed detector bank. The acces- 
sible energy range was Q— dependent. For the smallest 
Q energy transfers htu were simultaneously discerned be- 
tween — 1.5 meV and H-1.5 meV with a high resolution of 
hAoj < ±30 yueV. For the largest Q energy transfers hut 
were ascertained between —1.5 meV and -1-8 meV with 
a similar good resolution. Scattering data were taken at 
various temperatures T to elucidate its effect on the dy- 
namics. 

For BL scattering data were collected in a transmission 
geometry with the sample normal parallel to the inci- 
dent neutron beam. Stochastic motions were studied 
at four temperatures above the bulk freezing point of 
Tfb = 180 K: 250 K, 230 K, 210 K and 185 K. At each 
temperature data were collected for at least 4 h to guar- 
antee good signal to noise ratios. 

For CL two different scattering geometries were used 
to discern between radial and axial motions of hexane 
molecules in the cylindrical pores. In the first geometry 
the angle between pore axis (wafer normal) and incident 
beam was $ = 135°. Therefore the main component of 
the probed wave vector transfers pointed along the sym- 
metry axis of the pores, creating a high sensitivity to 
axial motions. In the second, complementary geometry 
the angle was $ — 45°. Probed wave vector transfers 
were nearly perpendicular to the pore axis. Therefore 
the scattering experiment elucidated mainly radial mo- 
tions of hexane molecules in the pores. Scattering data 
were taken in the liquid regime above the bulk freezing 
temperature at 250 K, 230 K, 210 K and 185 K as weh 
as below the bulk freezing point at 165 K. Measurements 
of liquid hexane at T = 165 K were possible because 
of the sizeable shift of the liquid-solid phase transition 
upon spatial confinement. [4 21 a shift which is for hex- 
ane not less than 20 K [43 in 6 nm wide silicon pores. 
Each individual measurement took again 4 h. The liquid- 
solid transition itself in confined hexane was probed by 
a sequence of 'quick-shot' measurements recorded while 
cooling the sample down from 165 K to 10 K. In steps of 
1 K data were sampled in each case for 5 min. 
Complementary scattering data were taken to allow and 
facilitate a quantitative data analysis. Background data 
were recorded for an empty aluminum cell and an evac- 
uated aluminum cell containing the porous silicon plate. 
These data sets were used to adjust the experimental 
results for BL and CL with regard to scattering contri- 
butions not originating from the hexane pore condensate. 
Additional scattering data were taken from solid hexane 
at T = 10 K. Here the complete incoherent scattering 
signal builts up in an elastic line due to disrupted molec- 
ular motions at low temperatures and obtained data can 
be used for resolution corrections. 
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IV. EXPERIMENTAL RESULTS 

This section is devoted to present a subset of the scat- 
tering data, which were recorded for BL's and CL's in 
the three-dimensional parameter space of wave vector 
transfer Qj energy transfer fiuj and temperature T. The 
presentation focuses for BL's and CL's on representa- 
tive scattering data recorded for a wave vector transfer 
of Q = 1.00 A^^. In the case of CL's the presentation 
is abridged and only scattering data sensitive to axial 
motions in the pores ($ = 135°) are discussed. Scatter- 
ing data recorded for different momentum transfers or 
recorded for different scattering geometries ($ = 45°) 
could be discussed equivalently. They do not show qual- 
itatively different characteristics as the ones discussed 
below. 

Fig. [2^ shows the incoherent quasi elastic scattering sig- 
nal S{Q,uj) recorded for bulk hexane at various tempera- 
tures T for a fixed momentum transfer of Q = 1-00 A~^. 
Neutron intensities are shown as function of energy trans- 
fers huj between —1.5 meV and -1-1.5 meV. Shown data 
were corrected for background scattering but not for en- 
ergy resolution R{Q,u]). Neutrons which gain energy 
in the scattering process contribute to the intensity at 
TiLo > meV, neutrons which lose energy contribute 
to the intensity at fiuj < mcV. Observed neutron in- 
tensities are always maximal for zero energy transfer 
{hu} = meV), that is elastic scattering. Intensities decay 
symetrically around the elastic line (Tiuj = mcV). Posi- 
tive and negative half width at half maximum (HWHM) 
are identical apart from the sign as expected according 
to the principle of detailed balance for energy transfers 
huj << k^T. Quantitatively the signals half width at 
half maximum decreases from 0.2 meV at T = 250 K to 
0.07 meV at T = 185 K. This reflects qualitatively the 
expected mobility loss of the hexan molecules while ap- 
proaching the freezing point. 

Fig. [^b exhibits similar to Fig. [5^ the incoherent quasi 
elastic scattering signal measured for hexane confined in 
the cavities of porous silicon. Data sets were again cor- 
rected for background contributions. A correction for 
energy resolution was not applied. The scattering signal 
exhibits the same qualitative behavior as discussed above 
for the BL. The signals half width at half maximum de- 
creases from 0.09 meV at T = 250 K to 0.035 mcV at 
T = 165 K indicating slowed dynamics close to the freez- 
ing point. Half widths for CL's are at given temperature 
about 50% smaller than respective half widths for BL's. 
Molecular motions in the confined system are therefore in 
average slower than in the bulk system. A more detailed 
and quantitative interpretation of this important obser- 
vation requires, however, a theoretical approach which 
accounts for different kinds of stochastic motions of hex- 
ane molecules in BL's and CL's and must carefully con- 
sider the dependent energy resolution R{Q,oj) of the 
experiment. 
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FIG. 2: (Color online) Incoherent inelastic scattering intensity 
S{Q,(jj) as function of energy transfer hu and temperature T 
for a fixed momentum transfer of Q = 1.00 A~^. Data are 
shown for BL (a) and for CL (b) with $ = 135°. 



V. NEUTRON DATA ANALYSIS: THEORY 

Hexane molecules scatter predominantly incoherently 
(due to the large incoherent cross-section of the pro- 
tons). Therefore, we analyze in the following our data 
with a model for the incoherent scattering cross section 
5ic(Q,a;) which takes into account the different kinds of 
stochastic motions that may cause incoherent inelastic 
scattering in our experiments. 

In the BL, long range translational diffusion as well as 
more localized motions, that are rotational reorientations 
of entire molecules or fast intra-molecular motions cause 
contributions to the inelastic scattering signal. Incoher- 
ent inelastic scattering due to translational diffusion is 
commonly described by a single Lorentzian term with 
g-dependent HWHM 7t(Q).141 It is: 



'(Q,^) 



1 



7t(Q) 



(1) 



For continuum diffusion as governed by Fick's second law, 
it can be derived that 7t = TiDQ^ ,where D refers to the 
so called self-diffusion constant. [4J] A more sophisti- 
cated, microscopic approach based on a jump diffusion 
model confirms the linear relationship between 7t and 
for small Q-values, but predicts deviations for larger 
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Localized molecular reorientations and intra-molecular 
motions can approximately be governed by [3 



1 



7loc 



Ao(Q)5H+(i-^(g))--^ 



(2) 



Ao{Q) is the so called elastic incoherent structure factor 
and emerges from spatial restrictions by definition im- 
posed on localized stochastic motions. [43] The HWHM 
7ioc of the Lorentzian term in Eqn. (5) might be related 
to a characteristic relaxation time r of localized motions 
by 7ioc = fi/T with T = Tge"^*/'"'^"^. Ea is the activa- 
tion energy of the considered processes. Assuming an 
independent superposition of long range diffusional pro- 
cesses and localized motions leads to the total incoherent 
scattering cross section: 



,A{Q) 7T 



(3) 



C* ( 
1 



MQ) 7T + 7loc 



■7^ 



Here (E) denotes a convolution in w and C is a Q- 
indcpendent scaling factor. 

In the CL the stochastic motions will be not qualitatively 
different, but at least two different molecular surround- 
ings have to be discerned. Molecules close to the pore 
walls are expected to exhibit different dynamic properties 
than molecules in the pore centre. The most obvious 
approach would therefore be to compose a total scatter- 
ing function as sum of two terms, each similar to Eqn. 
which represent wall layers respectively capillary conden- 
sate. Such an approach, however, would inevitably lead 
to a multitude of free and strongly correlated parame- 
ters in the model function. Consequently we propose a 
simplified approach to reduce the number of parameters: 



5r^'(Q,w) = C*((l-/c)( 
l-A{Q) 



AiQ) 



7t 



7t + 7ioc 



■7^ 



h uj"^ + (7t + 7ioc)^ ' 



(4) 



Here it is assumed that a fraction I — /c of the CL can 
be described similar to a BL as sum of two Lorentzians, 
while the complementary fraction /c of the CL is com- 
pletely immobile and is therefore represented by a delta 
function 5{ui) in the scattering cross section. 
A convolution in w— space of Eqn. ^ respectively 
Eqn. (|4]) with the experimental resolution function 
R{Q,oj) allows to relate theoretical models 



IiQ,uj) = R{Q,uj)®Sir\Q,uj) 
and observed scattering signals for BL and CL. 
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VI. NEUTRON DATA ANALYSIS: RESULTS 

Nonlinear least square fitting based on the discussed 
models was utilized to analyze the scattering data 
quantitatively in the probed T— and Q^range. Eqn. ([5]) 
and Eqn. ([5]) were adjusted to the BL data. In the 
fitting approach the HWHM 7t of the first Lorentzian 
in Eqn. ([3]), the characteristic relaxation time t for fast 
and localized motions, the elastic incoherent structure 
factor A{Q) as well as the scaling factor C were treated 
as free parameters. 

Fig. [3^ illustrates representatively the success of the ap- 
plied fitting approach for a particular set of Q and T. In 
the semi-logarithmic plot symbols exhibit the incoherent 
scattering signal S{Q,u}) taken at Q = 0.49 A^^ and 
T = 185 K corrected for background contributions. The 
solid line shows the best approximation of our model 
(Eqn.'s ([3])- ([5])) to the data. Different contributions 
to the model function ^ are shown as dashed lines. 
The Lorentzian term with HWHM convoluted with 
the resolution function is shown as short-dashed line. 
For small energy transfers {uj < 0.2 meV) this term 
is the main contribution to the scattering signal. The 
long-dashed line represents the second Lorentzian with 
HWHM 7t + 7ioc convoluted with the resolution func- 
tion. It is the main contribution to the scattering signal 
for large energy transfers uj > 0.5 meV. The resolution 
function itself is shown as a dot-dashed line. Fitting of 
scattering data at other points in the parameter space Q 
and T gives similar good results. A detailed discussion 
of the obtained fitting parameter and their dependence 
on Q and T will be given in the subsequent paragraphs. 
In the fitting approach the parameter C does not show 
a systematic dependence on the wave vector transfer 
Q. Independent of temperature T found values for C 
arc close to unity. Deviations from unity are smaller 
than 10% and as a direct consequence of Eqn. ([3]) and 
Eqn. ([5]) the in w— space integrated intensity of the 
scattering signal at given Q is considered to be equal 
to the integrated intensity of the respective energy 
resolution function R(Q,uj). The latter is in so far 
not surprising as we used the scattering signal from 
the solidified BL as resolution function in the fitting 
approach and not the scattering signal from standard 
vanadium samples. Constance of C illustrates only 
the conservation of incoherent scattering intensity at a 
particular Q. It is therefore a first successful consistency 
check for our model. 

Fig. 2^ shows the HWHM 7^ of the first Lorentzian in 
Eqn. ^ as function of the square of the wave vector 
transfer Q for probed temperatures between 250 K and 
185 K. For small wave vector transfers < Q/A~^ < 1 
the HWHM jt is at given temperature T proportional to 
as expected for translational continuum diffusion. [43| 
Linear fitting in this restricted Q— range was exploited 
to extract self diffusion coefficients Dt according to 
7t = TlDQ^. Obtained temperature dependent diffusion 
coefficients Z3t arc listed in the Fig. 2^ and Table HI The 
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FIG. 3: (Color online) (a) Symbols (circles) exhibit the in- 
coherent inelastic scattering from BL at T = 185 K. The 
solid line shows the result of nonlinear square fitting utiliz- 
ing Eqn. (131) and Eqn. The short dashed line represents 
the first Lorentzian term in Eqn. ^ convoluted with the en- 
ergy resolution function denoted as dot-dashed line. The long 
dashed line represents the second Lorentzian term in Eqn. Q 
convoluted with the energy resolution function, (b) Symbols 
(circles) exhibit the incoherent inelastic scattering signal from 
CL at r = 165 K ($ = 135°). First and second Lorentzian 
term of Eqn. @ were convoluted with the energy resolution 
function and are shown as short and long dashed lines. The 
dotted line illustrates the 5-function in Eqn. Q convoluted 
with the energy resolution. The solid line is the result of the 
fitting approach based on Eqn. Q and Eqn. (15}. 



diffusion coefficients for BL vary between 1x10 ^ cm^s ^ 
at 185 K and 3.1 x 10"^ cm^s"! at 250 K. Diffusion 
coefficients between 0.5 x 10^^ cm^s^^ at 185 K and 
2.15 X 10~^ cni^s"^ at 250 K tiave been found for bulk 
hexane in NMR spin-echo measurements performed by 
Douglass and McCall 1958 and 1959.(11 lH 
Deviations form the linear relationship between 7t and 
at larger scattering vectors Q arc expected. They 
mark length scales in/Q on which molecular diffusion 
can not be described any more by macroscopic models 
(Pick's law) but more microscopic approaches [i^ are 
required. 

In the fitting routine the paramater r was considered 
independent of Q as motivated by Dianoux et al..[50j 
Obtained r values suggest a time constant tq — 0.4 ps in 
the pico-second range and very low activation energies of 
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FIG. 4: (Color online) Half width at half maximum 7t of 
the first Lorentzian term in Eqn. ^ respectively Eqn. (Q 
as function of and T. Data are shown for isotropic BL 
(a), for CL with $ = 135° (b) and for CL with $ = 45° 
(c). Different symbols (black triangle, square, star, circle, 
grey triangle) refer to different temperatures (250 K, 230 K, 
210 K, 185 K, 165 K). Solid lines are linear fits to estimate 
the translational diffusion constant Dt. The fitting range was 
restricted to < / kr^ < 1. 



TABLE I: Translational diffusion constants for isotropic bulk 
hexane (Dbl) and for confined hexane along, respectively per- 
pendicular to the pore axis (Dcl , -Dcl) measured at different 
temperatures T. 



T [K] 


Dbl [cm^s~^] 


Dhf [cm^s-i] 


Dhl [cm^s-i] 




250 


3.1 X 10~^ 


2.1 X 10~^ 


1.8 X 10-^ 


0.68 


230 


2.4 X 10"^ 


1.7 X 10"^ 


1.6 X lO"'^ 


0.71 


210 


1.7 X 10"^ 


1.2 X 10"^ 


1.4 X 10"^ 


0.71 


185 


1.0 X 10"^ 


0.8 X 10'^ 


0.8 X 10"^ 


0.8 


165 


N.A. 


0.5 X 10"^ 


0.5 X 10"^ 


N.A. 



less than 10 meV for the localized motions governed by 
the second Lorentzian in the fitting function. Therefore 
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discussed, probably intra-molecular motions appear 
to be one to two orders in magnitude faster than 
translational diffusion with T— dependent relaxation 
times tt =< P > /6D between 20 ps and 70 ps. Here 
< P > denotes the mean square displacement of a 
hexane molecule in the time tt and was assumed to be 
given by the square of the nearest neighbor distance 
(c?NN = 0.65 nm) in the liquid. 

The data analysis did not reveal any significant tem- 
perature dependence of the elastic incoherent structure 
factor A{Q) at a fixed wave vector transfer Q. The 
obtained Q— dependence is shown in Fig. [S] (stars). 
A{Q) decreases continuously from 0.75 at Q = 0.5 
to 0.4 at Q = 1.5 A~^. The solid line exhibits an ex- 
pected decay in A{Q) for rotational diffusion of hexane 
molecules around their symmetry axis.fsoj Although 
it shows the same qualitative behavior as the data a 
quantitative agreement is not achieved. This is in so 
far not surprising as a multitude of different localized 
motions in the hexane molecule define A{Q) not only 
rotations around the symmetry axis. Further multiple 
scattering effects as outlined by Zorn et al. [HH might 
lead to this discrepancy between experiment and theory. 
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FIG. 5: (Color online) Elastic incoherent structure factor as 
measured for bulk hexane (stars) and confined hexane (cir- 
cles). Expected elastic incoherent structure factor (solid line) 
for rotation of hexane molecules around their symmetry axis. 

Eqn. (HI and Eqn. (O facilitated a quantitative analysis 
of the CL data. Free parameter in the fitting routine 
were the HWHM 7t, the HWHM 710c, the Elastic Inco- 
herent Scattering Factor A{Q), the scaling parameter C 
and the fraction fc of immobile hexane molecules in the 
pores. 

Fig. exhibits the excellent agreement found between 
CL data ($ 135°) and model function (Eqn. (g])) 
for a selected combination of wave vector transfer 
Q = 0.49 and temperature T = 165 K. Scattering 
signal S{Q,uj) and best approximation based on our 
model are shown, as well as the different contributions 
to the model function that were discussed above. As 
for bulk hexane translational diffusion and fast localized 
motions contribute significantly to the model function. 



But in order to facilitate a good approximation of the 
scattering signal at small energy transfers hui a per- 
centage /c > of hexane molecules is set in the fitting 
routine to be immobile. Scattering from this immobile 
molecules is shown as dotted line in the figure. It is a 
vital contribution to the model function in an energy 
range around the elastic line (huj = mcV) which 
corresponds to the half width half maximum of the 
energy resolution function R{Q,uj) at given Q. Applying 
our fitting approach to data that were recorded at other 
points in the parameter space Q, T, respectively for the 
second probed scattering geometry (<I> = 45°), leads to 
similar compelling results. A quantitative and systematic 
evaluation of the free parameters obtained by fitting the 
data will be given in the subsequent paragraphs. 
In the fitting routine parameter C exhibits the same 
quantitative behavior as already observed in the analysis 
of the BL data. Independent of wave vector transfer Q 
and temperature T values for C were found which do 
not deviate by more than 10% from unity. The elastic 
incoherent structure factor A{Q) found for the CL does 
not show any significant temperature dependence. Its 
Q— dependence in Fig. [S] (circles) exhibits exactly the 
same characteristics as observed for the BL. Paramter r 
predicts again relaxations times in the order of 1 ps and 
activation energies below 10 meV for the fast localized 
motions. 

Fig. |4}d and Fig. \^ exhibit the relation ship between 7t 
and the square of the wave vector transfer Q for various 
temperatures T. Frame b relates to scattering data 
sensitive to motions along the pore axis ($ = 135°). 
Frame c relates to scattering data sensitive to radial 
motions in the pore ($ = 45°). The linearity between 
7t and in the wave vector range < Q/A~^ < 1 was 
exploited to extract translational self diffusion constants 
Dt- Resulting diffusion coefficients for the different 
scattering geometries arc listed in the respective figures 
as well as in Table HI For $ = 135° diffusion coefficients 
from 2.1 X 10"^ cm^s^^ to 0.5 x 10"^ cm-^s^^ are found 
in the temperature range between 250 K and 165 K. For 
$ = 45° diffusion constants between 1.8 x 10~^ crn^s"^ 
and 0.5 x 10~^ crn^s"^ are found. 

Self diffusion coefficients in the mobile part of the 
pore condensate are about 20% to 30% smaller than 
respective ones in liquid bulk hexane. In the margin 
of error a directional anisotropy can not be ascer- 
tained. Translational diffusion along the pore axis is 
not quantitatively different from translational diffusion 
perpendicular to the pore axis. Radial confinement on 
a length scale of 6 nm (diameter of the pore), that is 9 
times the length of an hexane molecule in an all-trans 
configuration, seems not to be small enough compared 
to the approximated displacement V< P > « 0.65nm of 
a single hexane molecules in the time tt to affect radial 
diffusion coefficients . 

The Arrhenius plot in Fig. IH] exhibits self diffusion 
constants Dt reported in this and various other articles 
for bulk hexane and confined hexane as function of in- 
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FIG. 6: (Color online) Arrhenius plot for the T-dependence of 
the translational diffusion constant Dt' Squares denote bulk 
data. Triangles and stars refer to diffusion constants along the 
pore axis ($ = 135°) respectively perpendicular to the pore 
axis ($ = 45°). Solid lines illustrate linear fits to extract the 
activation energies Ea- Dashed lines exhibit IQENS results 
from Baumert et al. as shown in [s^ for mesoporous silica gels 
with 2 nm, 4 nm and 6 nm pore diameter and results form 
McCall et al. obtained in NMR spin-echo experiments. [49| 
Results of Baumert et al. are not corrected for a missing 
factor of 2n. 



Cooling 
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FIG. 7: (Color online) Temperature dependence of the 
amount fc of immobile hexane in the pores. Symbols rep- 
resent values extracted from the scattering data. Lines help 
to guide the eye. 



of 2 nm, 4 nm and 6 nm. The reported constants are 
significantly different from our results. However, this 
difference can be traced back to a missing factor of 27r 
in the calculations in the before mentioned article. The 
reported activation energies for tranlational diffusion 
in [36| are considerably larger than the one found 
here. This trend is partially explained by the pore-size 
dependence of the activation energy as illustrated in. (36j 
An analysis of scattering data in the temperature range 
between 250 K and 40 K revealed a temperature depen- 
dent fraction /c of hexane in the pores (Fig. [7]), which is 
immobile on the probed time scale of 1 ps to lOOps. At 
250 K about 5% of the CL showed no dynamics. This is 
significantly less material as required to entirely cover 
cylindrical pore walls with radius i? = 3 nm. Therefore 
it seems likely that this fraction of hexane represents 
molecules that arc trapped in deep adsorption sites of 
the less then perfectly smooth, but rough internal pore 
walls. With decreasing temperature fc increased at first 
slowly. Close to T = 160 K, that is the free zing point 



of hexane in 6nm silicon pores as reported in [43| about 
20% of the CL appeared to be immobile. This agrees 
fairly well with the amount of hexane on the pore walls 
at onset of capillary condensation and makes a strong 
case that hexane molecules in close vicinity of the walls 
lose their mobility while the CL approaches the freezing 
point. This population corresponds likely to a flat- lying, 
strongly adsorbed monolayer of hexane molecules, which 
sticks to the pore wall, similarly as it has been inferred 
from ellipsometric and rheological measurements for 
medium-length n-alkanes in the proximity of silicon and 
silica surfaces, [l^ HO, HE Here it should however 
be noted that between 180 K and 160 K hexane already 
solidifies in a small number of pores with a larger than 
average diameter (> 6 nm). These immobile hexane 
molecules will also contribute to fc as the scattering 
experiment can not discern between hexanc stuck on 
the pore walls and hexane solidified in larger than 
average pores. Below 160 K the immobile fraction 
increases rapidly. The hexane solidifies in the pores. 
This solidification process spreads over an temperature 
range of 20 K associated with the pore size distribution. 
Below 130 K there is no mobile hexane left. The fraction 
fc of immobile hexane was extracted from data which 
were recorded while the sample was cooled down. 



VII. CONCLUSIONS 



verse temperature 1/T. Symbols illustrate the diffusion 
constants reported in this article for BL and the mobile 
part of the CL. Activation energies as estimated from the 
arrhenius plot are similar for BL and CL. The obtained 
value of 60 meV— 70 meV is 20% smaller than the one 
reported by McCall. [i^ Fig. [6] shows also Diffusion 
constants reported by Baumert et al. [36| for hexane 
confined in mesoporous silica gels with pore diameters 



We performed timc-of-flight experiments at spectrom- 
eter IN5 (Grenoble) to probe stochastic motions in liquid 
bulk hexane and liquid hexane confined in porous silicon 
which exhibits almost cylindrical pores parallel-aligned 
along the surface normal with an average pore diameter 
of 6 nm. 

Long range as well as localized motions of hexane 
molecules confined in the porous silicon matrix were 
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found to be slightly slower than in the BL. A temper- 
ature dependent fraction /c of the CL was found to be 
immobile on the probed time scales. At 250 K about 
5% of the confined hexane, did not show any dynamics 
at all. These 5% of hexane likely accounts for molecules 
pinned in deep adsorption sites of the rough pore walls. 
At 165 K before onset of freezing about 20% of the pore 
filling, that is the amount of hexane on the pore walls 
before onset of capillary condensation, was immobile. 
In the mobile part (1 — /c) of the CL translational 
difi^usion constants Z^t were reduced by 20% to 30% 
compared to the BL, but not by an order of a magnitude 
as inadvertently claimed by Baumert et al..(36j A reason 
for the slowed dynamics in the mobile hexane can not 
entirely be elucidated by the performed experiments 
and analysis. It might be a direct result of spatial 
confinement imposed on molecules in the pores or to 
some extent an artifact of the simplified fitting approach 
which only discerns immobile hexane molecules and 
mobile hexane molecules. A radial gradient in the 
diffusion constants, that is bulk diffusivity in the pore 
centre and a decreasing diffusivity closer to the pore 
walls, seems to be intuitive and suggested both by 
Molecular Dynamics simulations and experiments. 0, [53| 
But attempts to apply models that account more clearly 
for a radial dependence of dynamics failed due to a 
multitude of strongly correlated fitting parameters and 
were therefore not suited to perform a reliable and 
meaningful analysis of the scattering data. 
Follow-up experiments should exploit partial fillings of 
porous hosts with hexane to probe radial variations in 
diffusivity directly. Previous such studies employing 
mesoporous MCM-41 silica matrices indicated, however, 
that additional to the stochastic motions of the molecules 
in the pore centre new diffusion mechanisms, i.e. at the 
solid/vapor and at the liquid/vapor interfaces of the 



partially filled matrices, have to be considered. (33|, |5^ 
They may hamper an extraction of the radial depen- 
dency by such filling-fraction-dependent measurements. 
Nevertheless, probing CL's before and after the onset 
of capillary condensation in the pores should allow to 
discern more clearly the dynamics of hexane molecules 
close to the pore walls and dynamics of hexane molecules 
close to the pore centre. 

An anisotropy in diffusion coefficients was not found 
in the CL. Translational diffusion of hexane molecules 
along the pore axis was not different from diffusion 
perpendicular to the pore axis. A radial confinement 
on 6nm was not sufficient to affect the diffusion of 
hexane molecules which is characterized by a mean 
displacement of about I = 0.65 nm on a time scale 
of tt- This finding is corroborated by high-resolution 
optical birefringence measurements on hexane and 
other rod-like molecules of comparable size confined 
in cylindrical nanochannels, whic h g ave no hints of 
any static orientational order, (ssi ISo l Experiments on 
confined liquids in porous hosts with less wide pores are 
expected to show a different outcome. If the diameter 
of the pores gets closer to the average displacement I 
of hexane molecules in the time tt the diffusion in a 
cylindrical pore is expected to become more and more 
one dimensional and pref erentially directed along the 

pore axis. H [13, [H-dl [53 
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